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> ■ Abstract 

We show that scalar, 0-form, Galileon actions — models whose field equations contain only second 
C\) ■ derivatives — can be generalized to arbitrary even p-forms. More generally, they need not even 

depend on a single form, but may involve mixed p combinations, including equal p multiplets, where 
odd p-fields are also permitted: We construct, for given dimension D, general actions depending 
(^i on scalars, vectors and higher p-form field strengths, whose field equations are of exactly second 

derivative order. We also discuss and illustrate their curved-space generalizations, especially the 
delicate non-minimal couplings required to maintain this order. Concrete examples of pure and 
mixed actions, field equations and their curved space extensions are presented. 
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I. INTRODUCTION 

The geometric ancestors of Galileons llH5| are the Gauss-Bonnet-Lovelock (GBL) actions 
/= j d D xe^-e^-R^ aP ...R..e pl ...e..., (1) 



powers of the curvature R whose field equations are nevertheless independent of higher than 
second metric derivatives. This is achieved by virtue of the Bianchi identities, due to which 
the /^-variations do not contribute; only the explicit vielbeins' do, as is especially clear in 
vielbein/spin connection formalism. Here R is the "field strength" of the (non-Abelian) spin 
connection o; MQ ,^(e); the Levi-Civita symbol e^— is a tensor density, while e 01 "- is a world 
scalar; (//, u, . . .) and (a, (3, . . .) are world and local Lorentz indices respectively. These 
actions are dimension-dependent, yielding vanishing field equations below a certain D, such 
as D = 5 for R 2 and D = 7 for R 3 . More explicitly, for D = 5 say, one e Ma is required to 
contract the two leftover indices in {eeRR)^ a 1 while there is no e^ a , hence no field equation, 
in D = 4. The mechanism is simple, and as we shall see below, universal: First note that 
8Rfiua/3 = T)hi8wi/]aP) where V is the usual covariant derivative with respect to the spin 
connection (acting also on local indices), and Sco is a world vector. Therefore, integrating T> 
by parts (freely past all vielbeins of course) onto the remaining Riemann tensor(s) gives by 
the cyclic Bianchi identities. So the GBL field equations, (eeR . . . _R) Ma = 0, just result from 
removing (any) e^ a in ([1]) and are manifestly independent of higher than second vielbein 
derivatives. 

Galileons are scalars whose field equations depend only on second derivatives, hence are 
invariant under constant shifts of the fields ("positions") and their gradients ("velocities"), 
recalling old-fashioned Galilean invariance. [Note that this invariance is only meaningful 
in flat space, since there are no constant vectors or tensors in curved backgrounds.] Their 
actions bear a formal resemblance to the GBL systems, when expressed as jj] 

1= f d D xe^-£ a P-d fl <Kd a Tr{d v dp<ir)...{ddTr). (2) 

Again, the variations only leave second order equations (ddir) . . . (ddn) = 0, and for suffi- 
ciently low D, where the dirdir ~ "e" are absent (for a given power of ddn ~ "-R"), Eq. (T5]) 
has vanishing variation. This (slightly imperfect) similarity led us to conjecture that ([2]) 
could be obtained from a GBL-like action usinga metric suitably parametrized in terms of 
dn; it was indeed elegantly confirmed recently [5| for the R 2 case (and in a suitable limit), 
as a byproduct of a brane analysis. 

The purpose of this Letter is to generalize the above models by noting that the properties 
of 0-forms underlying (J2J) are actually shared by arbitrary even p-forms, and extend to 
any (dimensionally allowed) admixtures of various p-level fields. Surprisingly, we found a 
fundamental divide between even (scalars, . . . ) and odd (vectors, . . . ) models. The latter 
turn out, despite initial appearances, to have empty flat space actions 1 (except of course 
the standard Maxwell-like C = F 2 ), i.e., devoid of field equations for any (D, p — 2n + 1). 
However, as we discuss below, they may appear in mixed form, or in multiplets of single 
p-form, models. 



1 Nevertheless, covariantized versions of trivial flat space actions can produce nontrivial field equations, 
proportional to the curvature, as will be illustrated in Sec. IIVI 



We will work primarily in flat space in order to focus on our main results. As for scalars, 
the key ingredient here is that the forms'"field strengths" uj p+ \ = dA p , are curls which do 
not become covariant; only the explicit V in Vw does. Using these "gauge-invariant" field 
strengths rather than ordinary gradients is both essential to the Galileon aspect and excludes 
their, possible ghost, lower spin gauge components. 

Retaining second order upon extension to curved backgrounds is nontrivial; even for 
scalars, the minimal coupling extension of fl2]) gave rise to third derivative terms in its 
stress tensor and hence in the associated gravitational field equations, as well as to third 
metric derivative terms oc V-R in the n equations. A delicate set of additional, non-minimal, 
couplings, involving the full curvature tensor was required [3j, |4| to remove these. This 
program, though correspondingly more complicated, can in fact be carried out for our present 
generalized framework in a fashion similar to that of |4| for scalars. Instead of detailing here 
the straightforward but still rather lengthy derivation for the most general case, we will 
display some examples of successful covariantization in Sec. IIV1 



II. p-FORM ACTIONS 

We start, to emphasize the pitfalls in this problem, with the obviously simplest — but 
actually empty — generalization from a scalar (j2j) to a one-form A p with field strength 

F^u = d[^A v y. 

1 = J d D xe^-e a ^-F, u F a ^d p F l5 . . .)(d e F UT ...), (3) 

where the parentheses contain products of OF and indices are connected as follows: In 
the first parenthesis, the index of the derivative d is contracted with the first e^— 
whereas those of F are contracted with the second e a/3 '", and inversely in the sec- 
ond parenthesis. The integrand of ([2D is a total divergence which we may write as 
l d ^ £ nv... vx ... £ ap...ec... F^F aT F a p(d p F lS • • • )(d c F vx ...)]. This equality follows by noting that 
|<9 e manifestly annihilates all its operands but F pv F ar (on each of which it acts identically), 
where its actions reproduce the original Lagrangian. Since this conclusion is due to the 
evenness of F^F^ upon exchange of their indices \iv -B- ar, the difficulty obviously persists 
for all (D, p — 2n + 1). However, as we shall see in the next section, odd p models can be 
revived if they are allowed to depend on more than one A^. 

Fortunately, the direct extension of ([2]) does exist for even p; it formally resembles ([3]) 
where now u)\ pu ,,, = d[\A^ u ^]. In detail, the general (now non- vanishing) action is 

d D xe^-e ap -u^...uj a p..{d p uJ 1 s... ■ ■ ■ )(d e u aT ... ...). (4) 

As in Eq. ([3]), when the derivative d of a gradient du is contracted with one of the two e, the 
indices of u p+ i must be contracted with the other e, otherwise the action would vanish by 
virtue of the Bianchi identities (i.e., [d, d] = 0). The two parentheses of (jlj) must contain the 
same number of terms, not greater than (D — p — l)/(p + 2). [In the lowest, p = 0, case, an 
odd total number of du is also permitted, cf. ([2]).] For fewer terms, there remain extra open 
indices on each e that must be contracted between them (using vielbeins in curved space, of 
course); this yields (up to an overall factor) the same Lagrangian as in the lowest possible 
dimension. There is actually a formal resemblance between (jlj) and ([2]) that emerges from 
considering the (p + 1) indices of the form u as a multi-index M replacing the single index 



of 7r M . Then, just as each of the two indices of d a d,/K must be contracted with different e, 
here a and M of d a UM must belong to different e. 

It should be clear from our notation and the Bianchi identities that the field equations 
depend homogeneously on du and not at all on u, hence they enjoy the corresponding 
Galilean invariance, this time under a shift of (A pu , u^ vp __) by constants antisymmetric 
tensors (c^.i, fc[ w ...]). For completeness, let us run through the argument, entirely akin 
to those for gravity and scalars: First, if either pure u is varied, the explicit curl on its 
A can only land on the other u, since Bianchi annihilates any d p d[ a U/3^,,,] (or <9 a <9[ M u;„p...]). 
Likewise, varying any du ~ d p d a A factor forces each of those two d to land on one of the 
two pure u; all other landings vanish, again by Bianchi. The simplest version of Eq. fll]), 
valid for D > p + 1, does not contain any derivative du, and is thus the standard kinetic 
term u 2 (valid for all p of course, though dynamically non-trivial only when D > p + 1). 
The first novel p-form Galileon action, involving just two du factors, requires D > 2p + 3. 
For p — 2, it reads explicitly 



d 7 x e^vx^secn ^ ^ g^^ q^ 



T<PX 



36 / d 7 x 



-9(u^ a u^u T ^ x u^) - \%{u\"u p Ju^^u^ p ) 

-3(co^ x co p(TT , x ) 2 - 9(cu^co p(TT , x ) 2 + 18(uvX CT V) 2 + 9(w^a;^ r ) 2 
-%u pvp ^\ J 2 - (u^uv^) 2 + (cu 2 )(cu^ a ) 2 - 3(co 2 )(co^ !p ) 2 ] . (5) 

Its field equation 

^T^aMSetv dpUah Q aU&< drtUrm = (6) 

is obviously of pure second order; we do not display its, 23 term, expansion. 

We conclude this section by discussing another amusing (if somewhat tangential) parallel 
between tensors and forms, which evokes the well-known conversion [6] of pure divergence 
D = 4 GB into general relativity (GR), with or without cosmological term: upon adding 
±A.e pa e v p to each R^ap in the topological invariant (PQ), it becomes proportional to the GR 
action, since the cross-term in ee(R± Aee) 2 is the scalar curvature, while the A 2 term is the 
volume density, cosmological, term. Subtracting (i? + Aee) 2 and (R — Aee) 2 actions removes 
the latter. For scalars, we similarly add im 2 ^^ to each of the two d^daTc in the pure GB- 
like / = J eeddir ddTt: this leads to the (massive or massless by subtraction) Klein-Gordon 
action. These extensions can be made for all forms: thus for the vector Proca/Maxwell 
actions, add ~ im 2 ?]^^^] to each d p F a p in the otherwise vacuous action I = J eedFdF, 
etc. Note that our mass construction is valid in all D [> (p + 2) of course]; that of GR, for 
all D > 4. That the former has a curved space extension is also obvious. 

III. MIXED FORM ACTIONS 

Our actions (jlj) can be further generalized by including several species, i.e., mixtures of 
various unequal p-forms compatible with a desired D. Labelling these species by (a, b, . . . ), 
the action takes the formal expression 

I = jd D xe^-e^u; u ^{d P ^s... ■ ■ ■ )(^ T ... . . . ). (7) 



The number of indices contracted with the first and second e must be the same and not 
greater than D, but the two parentheses may now involve different species and therefore 
a different number of terms. Here, Bianchi again ensures (exactly as for the single species 
version) that only duo appears in the field equations. Hence (flat space) Galilean invariance 
under translation of all (A p , co p+ i) by constant antisymmetric tensors (c p , k p+ i) is preserved. 
Note that odd forms are also allowed in (j7j), subject to various symmetry constraints. For 
example, no more than two OF factors can be present, otherwise (J7J) would involve at 
least one product of the form d^F^dpFjs, where the two d and the two F are respectively 
contracted with the same e tensors. Hence their indices can be interchanged by three 
permutations, /x -B- u, a(3 -H- jS, so they vanish identically. This single dFdF ceiling 
obviously also applies to higher odd p-forms; instead, the p = 2n models, being even under 
such permutations, may contain arbitrary powers of du consistent with a given D (but 
conversely, see below for limitations on even p actions). 

Let us quote two simple, mixed & 1-form, nontrivial examples; the first Lagrangian is 
defined in any D > 3: 

C = e^ p e a ^ F^Frf d p d^ = 4F"*F> i/a , - 2F 2 Dvr. (8) 

Both its 7r and A\ field equations are obviously of pure second order; explicitly, 

(F mp f - 2(i^J 2 = 0, (9) 

F^^tt^ + F^tt'^-F^Dtt = 0. (10) 

Similarly, in D > 4, the mixed model 

c = £ ^e a ^ s d^d a nd u F^dsF pa 

= -8(7r^F^F p ^7c^) + 4(tt^F wp ) 2 + 2{ir>"F vp>(t ) 2 

-4(vr,,F^ J 2 - 2{^ p f{F vp ,f + 4(7r, M ) 2 (F^) 2 (11) 

also yields pure second order -n and A\ field equations: 

±{>K, pv F»^F po *) - 2(F^7r^F^ a ) + 2(F^ p rr^F^ a ) 

-(F pa ^F<>%) + (Dn)(F^ p ) 2 - 2(D7T)(F^J 2 = 0, (12) 

2(tt,, p F V P ) + 2{ir»F lu „ir">) + 2(^\ p ,F^ >v ) - (n 4W ) 2 (F\) 

-2(Ok)(^F x ^) - 2(OK)(Tr\F^ u ) + (Dvr) 2 ^) = 0. (13) 

An even simpler class of mixed actions involves a single p-order species, but now as a 
"multiplet" A a , for instance pure scalars but with different 7r a replacing the single one 
in 02]). This extension even resuscitates odd-p actions: For instance, the simplest bi- vector 
Lagrangian of the type ©, £ = e ^P UT £ a ^ 5e F pu F^p d p FKd e F^ TJ is obviously no longer a 
total divergence. Our reasoning below Eq. flU, showing that the field equations do not 
involve higher order derivatives, may also be generalized to non-Abelian gauge bosons A a 
and their field strengths F = dA + A A A, although both the invariances under constant 
shifts, A a p ->• A p + c p and F pv ->■ F pv + k? ,, would then be lost. Indeed, if V denotes the 
covariant derivative with respect to the internal space (like the T> below f$B) with respect 
to tangent space), then the Bianchi identities D, F° ■, = still hold, therefore Lagrangians 
of the form C = e^-e 01 ^- F^ v F^{V p F^ s . . .){V e F% T ...) define non-linear extensions of 
Yang-Mills theory, while keeping field equations of second (and lower) order. 



It is worth noting that one may also add undifferentiated powers of u beyond the two 
in the generalized models (171), provided all indices of any one u (whatever its p-order) are 
contracted with those of a single e tensor, but not "across" both. Also, no more than two 
undifferentiated even p-field strengths u°; +1 are allowed for the same species a, otherwise the 
action would vanish by oddness, while any number of odd p field strengths may be present. 
The same reasoning as above indeed shows that no higher derivative than du is generated 
in the field equations, i.e., that they depend at most on second derivatives of the p- forms 
A. On the other hand, these field equations now involve some pure u in addition to the 
usual du factors, because at most two derivatives are generated by varying the du terms 
of the action, so that they can act on at most two of the undifferentiated u. Therefore, 
this generalization with more than two pure u results in a loss of the "velocity" invariance 
u — > u + k. A simple D > 4 example of this type is 



^F^TT^TV) + 8(7^^7^,,) + 2(F 2 )(7r, M 7T 



/.iw 



7T, 



+A(ir^) 2 {F» a F p °<K. vp ) - 4(7r iM F^iVrO(D7r) - 2(7r iM ) 2 (F 2 )(D7r). (14) 

As is clearest from the first expression, its variations involve both first and second (but no 
higher) derivatives of it and A^. 

IV. GRAVITATIONAL COUPLING 

As stated in the Introduction, second-order preserving extension of even the scalar flat 
space actions to curved backgrounds was a rather complicated process, one that becomes 
more combinatorially involved for higher forms. We content ourselves here with giving 
the explicit non- minimal extensions for four of our cases ([S]), dHD, (TTTT) . (TT41) . that avoid 
higher derivatives in both the matter and gravitational (that is, through T^) field equations. 
These terms are constructed as for scalars in [4]: All possible pairs of gradients, du a du b , 
must be replaced by suitable contractions of the undifferentiated u a u b with the Riemann 
tensor, and added to the minimally covariantized flat-space action with suitable coefficients; 
somewhat more involved counting shows that they require factors oc (p a + l)(pb + 1), where 
p a fi denote the orders of the forms A a,b . One other difference in the p > construction 
is that V^u;^... are to be distinguished from V a u pu ..., essentially because of their different 
e-index contractions, a distinction irrelevant to the original scalar, 7r ;MQ , = 7r ;Q , M , case. One 
common feature is that flat-space Galilean invariance is also not restorable by consistent 
covariantization (nor should it be expected, absent constant vectors or tensors in curved 
space): the equations now necessarily depend on both second and first derivatives of the 
fields. For (jSJ), the added terms are: 

A/ = ~ j ' d 7 xe^ paT ^e a ^ 5 ^ u, up u aM u XaT u\ e R vxC , n 

= 54 f d 7 xy^[24(u^u XTp u TVX u^R^ pa ) + l2{u Xr ^u mp u^ R PP U(T ) 

+4(u 2 )(u x , u u x ^R p % a ) + 12(u^u^u r ^u TXa R;) + 18(u, UT u^u^ T u vxa R;) 
-\Q{u 2 ){u pup u, ua R a p ) - 3{u, vp u pu °) 2 R + (u 2 ) 2 r\ . (15) 



Similarly, the mixed D > 4 example (TTTj) acquires the terms 
AC = e^E a ^ 5 d^d a <KF Xv F\R pa „ l5 



-9 



-2{n^, u F pX F a x R^°) + 4(ir^F up R»°K ff ) + 2(^) 2 (F ua F/R^) 



+2(F 2 )(n, p R^n^ + {n^) 2 R - (k,) 2 (F 2 )r\ . (16) 

In contrast to the above models, the mixed D > 3 example (IE]) and the "non-Galileon" La- 
grangian f|T4|) actually require no additional terms to preserve second order, since they only 
contain a single vulnerable — because second order — ddir factor. It is clear by inspection of 
OH]) and (fi~4]) that all (covariant) third derivatives arising from variations here always have 
the form of a commutator [V, V] acting on a dA, that is a — harmless — curvature times 
first derivatives of fields. 

Our final model illustrates the observation made in our footnote that actions trivial in 
flat space can have non-trivial, dynamical, curvature-dependent extensions: Consider the 
vector models Q, or more generally actions (]!]) for any odd p, which are vacuous in flat 
space. Their minimal covariantizations are clearly both nonvanishing and of third order. 
However, one may also add appropriate non-minimal terms that both remove the offending 
higher derivatives and remain non-trivial. Indeed, the simplest case is the lowest Galileon 
D = 5 vector action, 

/ = ! d 5 xe^ p ° T e a ^F^F aP V p F l5 V e F aT 

= -i j 'd 5 xe^ p,rT e a ^ Se F^F^F^R^. (17) 

The last equality in ffT7]) exhibits the model's curvature-dependence, and is obtained from 
the first expression by parts integration. [The metric variation of the curvature in the second 
expression ffTT]) yields a non- vanishing T pu = d a dpH^ a ^ u ^ even in flat space, despite the 
model's triviality there; no paradox ensues since this pure superpotential form has vanishing 
Lorentz generators.] The third derivatives in the resulting field equations can be removed 
by adding the counter-term 

A/ = J d 5 x e^ T e a ^ 9e F pu F a(3 F x p F Xj R arSe . (18) 

It differs from the action ffTT]) itself simply by an overall factor and the index change 5 -H- A 
in the last two terms. Their sum, 

I + AI = -8 f d 5 xy/=g 4(F^F^F AT F rCT C M1/p(r ) + ^{F^F^F^F^C^ 

+{F 2 ){F^F^C fiVpff )\, (19) 

depends only on the Weyl tensor C^ upa (for no obvious D = 5 reason, though ( fl9l is 
manifestly conformal invariant in D = 10); as per design, both its T pv and field equations 
depend on at most second derivatives. 

Details of our models' constructions, of their general non-minimal compensating gravi- 
tational extensions, applications for instance in the spirit of [7], and other open questions, 
e.g., possible supersymmetrization, may be presented elsewhere. 
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